Hybrid Stars in an SU(3) Parity Doublet Model 
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We apply an extended version of the SU(3) parity model, containing quark degrees of freedom, to 
study neutron stars. The model successfully reproduces the main thermodynamic features of QCD 
which allows us to describe the composition of dense matter. Chiral symmetry restoration is realized 
inside the star and the chiral partners of the baryons appear, their masses becoming degenerate. 
Furthermore, quark degrees of freedom appear in a transition to a deconfined state. Performing an 
investigation of the macroscopic properties of neutron stars, we show that observational constraints, 
like mass and thermal evolution, are satisfied and new predictions can be made. 
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I. INTRODUCTION 

The study of strong interaction physics under extreme 
conditions is a central topic of nuclear physics with a 
large number of experimental and theoretical programs 
focusing on the area. These conditions comprise large 
temperatures, densities, as well as extreme values of nu- 
clear isospin. In ultra-relativistic heavy-ion collisions a 
very hot fireball is created in the collision zone and, at 
high temperatures, hadronic matter is assumed to melt 
into its constituents, quarks and gluons. The net baryon 
density in such reactions is determined by the beam en- 
ergy. At the planned energies in the upcoming FAIR 
(Facility for Antiproton and Ion Research) at GSI, a hot 
and relatively dense system will be produced. A central 
point of these investigations is the understanding of the 
quark-hadron phase transition. However, while relatively 
high densities and finite temperatures might be reached 
in laboratory experiments, the study of neutron stars is 
essential if one wants to probe the low temperature and 
high density region of the phase diagram of strongly in- 
teracting matter. 

Observations of neutron star masses (and possibly 
radii) provide the standard way of constraining the inner 
composition of these objects. More precisely, recent ob- 
servations have set new high mass constraints for neutron 
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stars (PSR J1614-2230 Q being the most important), 
and equations of state aimed to describe compact stars 
are thus expected to provide objects with high masses 
In the case of hybrid stars (neutron stars with a 
quark core in its center), a quark phase based on a sim- 
ple non-interacting quark model like the MIT bag model 
tends to reduce the maximum mass significantly (see the 
discussion in 0). A quark phase that includes strong 
repulsive interactions, however, may have an equation of 
state quite similar to a purely hadronic one. This pre- 
vents the softening of the matter and the drop in maxi- 
mum mass It has been found, however, that mod- 
els with a strong repulsive quark-quark interaction make 
the description of lattice results at /is = nearly im- 
possible [Ill4l3l |. since including a repulsive mean field 
interaction strongly decreases the quark number suscep- 
tibility. Therefore, most current successful hybrid star 
models would fail in this regard, if they were to be ap- 
plied to the high-temperature, low-density regime. 

At high temperatures, QCD exhibits a crossover to 
a deconfined phase and the quarks and gluons become 
the dominant degrees of freedom fl4| . The temperature 
at which this transition takes place is estimated to be 
T dec » 150- 160 MeV (IE El. A phase transition is also 
expected to take place at high densities, where baryons 
start to overlap, and quarks and gluons become the ef- 
fective degrees of freedom. This indicates that at some 
point a hadronic model will not be able to appropriately 
describe the matter present inside a neutron star, and a 
deconfinement mechanism needs to be introduced. 

In the following, we will discuss a theoretical approach 
that is able to describe the conditions found in compact 
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stars as well as those created in heavy-ion collisions. The 
aim of this approach is to find a unified description for 
the thermodynamic properties of QCD which can be ap- 
plied to compact stars and heavy ion collisions at differ- 
ent beam energies while being in accordance with lattice 
QCD results at vanishing baryon density. An approach of 
this type is essential if one wants to investigate a phase 
diagram that has a region of a cross-over transition to 
quark matter (as it is clearly established by lattice QCD 
calculations) and a first-order transition at high densities 
and low temperatures. This is not possible by combining 
separate models for the hadronic and quark phase. Such 
a first-order transition at low temperature has also been 
observed in recent lattice monte-carlo calculations of an 
effective QCD Lagrangian (T3] ■ In addition, a model that 
can cover the physics at high temperature as well as at 
high density can serve as an important tool for the in- 
creasingly important studies of black hole formation and 
black-hole neutron-star mergers, where temperatures up 
to 90 MeV and high densities might be reached [HI, [l9| . 

In a previous paper, an extended quark SU(3)f par- 
ity doublet model was introduced for this purpouse [2(| ■ 
In the parity model jU [22j , the mass splitting between 
the nucleons and the respective chiral partners is gen- 
erated by the spontaneous breaking of chiral symmetry 
and their coupling to the corresponding order parameter, 
the scalar (sigma) field. The same applies for the baryon 
octet in the SU(3)f case [23|, [24J. When chiral symme- 
try is restored and the sigma field vanishes, the chiral 
partner masses become degenerate. A good description 
of nuclear saturation properties, as well as neutron star 
observables can easily be achieved within this formalism 
[25l I2H ] . In Ref . [I?} a phase diagram for chiral symme- 
try restoration was calculated using an SU(2)f version of 
the parity model. 

In addition, the parity model has been shown to de- 
scribe the lattice results at \xb = and can be used 
in dynamical models for relativistic heavy ion collisions 
[28l - l30| . We will apply this model to the study of neutron 
stars. A complete phase diagram for iso-spin symmetric 
matter, using the extended version of the SU(3)f par- 
ity model, which contains quark degrees of freedom, has 
been calculated in [2(j. As expected, at low tempera- 
tures the nuclear matter liquid-gas phase transition and 
the chiral symmetry restoration are of first order. Within 
this approach the deconfinement phase transition, on the 
other hand, is a crossover. Furthermore, chiral symme- 
try restoration and deconfinement do not coincide, and 
at intermediate densities the matter is chirally symmet- 
ric but still confined. Whether such a chirally symmet- 
ric hadronic phase can be the N c = 3 equivalent of the 
N c = 00 quarkyonic phase [3l[ is still subject to debate 

In this paper we investigate the properties of 
electrically-neutral chemically-equilibrated matter in the 
framework of the model described above. We will inves- 
tigate the influence of the chiral partners, hyperons and 
quark matter on the macroscopic properties of a neutron 



star, such as its gravitational mass, radius, and thermal 
evolution. 



II. THE SU(3) PARITY DOUBLET MODEL 

Considering that the baryons are grouped in doublets 
{B + and B-), in which they belong to the same multi- 
plet, the components of the fields ip± and 1^2 transform 
under L and R rotations like 

ip[ R = Rip 1R , i>' 1L = Ltp 1L , 

1p' 2R = Llp 2 R, lp' 2 L = R ^2L ■ (1) 

This formalism allows the presence of a bare mass term 
in the Lagrangian density that does not break chiral sym- 
metry 

^0(^275^1 - ^175^2) = 

m (^2LlplR - ^2RlplL ~ 4>lLlp2R + $1R^2l) , (2) 

where too represents a mass parameter. Since the term 
proportional to too mixes the upper and lower compo- 
nents of the parity doublets, one diagonalizes the matrix 
by introducing new fields B with a diagonal mass ma- 
trix. Keeping only the diagonal meson contributions, the 
scalar and vector condensates in the mean field approxi- 
mation, the resulting Lagrangian reads 

i i 

+ ^{B i ^{g ui ^+g pi p^+g (t>i ^)B i ) , (3) 

i 

where the coupling constants for the baryons with the 
mesons to, p and the strange meson <p are gsu>, gBp, and 
gBcf>, respectively. We do not include a vector-meson self- 
interaction to avoid an extra softening of the equation 
of state, which in turn would yield neutron stars with 
very low maximum mass [26l [35| . This happens because 
as the coupling constant increases, the respective vector- 
isoscalar field decreases (in order to reproduce saturation 
properties) , thus making the repulsive part of the strong 
force less pronounced. 

The expression for the effective masses of the baryons 
for isospin symmetric matter reads 

±g%a±g i $Q. (4) 

In the above equation, the term that does not correspond 
to the bare mass is generated by the scalar mesons er, 
the strange £, and by the SU(3) breaking mass term 
with m s = 150 MeV. The last term is responsible for 
the generation of an explicit mass corresponding to the 
strangeness n s of the baryon. The values of the various 
coupling constants of the baryons with the meson fields 
are given in Ref. [2(|. The couplings g^ are chosen in or- 
der to reproduce the splitting of the masses of the parity 



partners in vacuum, further assuming, for simplicity, an 
equal splitting of the masses of the various baryonic dou- 

(2) (2) 

blets (g ai = —0.850, g^ = 0). The coupling constants 

are chosen in order to reproduce vacuum masses of 
baryons. The bare mass parameter is set to too = 810 
McV. Such a high value is necessary in order to reproduce 
reasonably massive stars while keeping the compressibil- 
ity of saturated matter at reasonable values [25|, [26| . The 
value for the strange quark mass m s is fixed to 150 MeV. 

The nucleonic vector interactions are tuned to repro- 
duce reasonable values for the nuclear ground state prop- 
erties (binding energy per baryon of ~ —16 MeV and 
baryonic density of po — 0.15 fm 3 at saturation) while 
the hyperonic vector interactions are tuned to generate 
reasonable optical potentials for the hypcrons in ground 
state nuclear matter (J7 A ~ -28 MeV, Uy, ~ -14.5 
MeV and U 3 ~ -18 MeV). The equations of motion 
are obtained by minimizing the grand canonical poten- 
tial (a function of baryonic chemical potential and tem- 
perature), and arc solved sclf-consistently in mean field 
approximation. For non-zero temperature calculations, a 
hadronic heat bath is included. 

We choose the vacuum mass of the chiral parter for 
the nucleons to be 1535 MeV, the lighter resonance with 
spin 1/2 and negative parity. Lower values for this mass 
(~ 1200 MeV for instance) would lead to chiral restora- 
tion takingplace at lower densities [35| , and to less mas- 
sive stars [25[ . In addition, using a formula for the width 
as a function of the bare mass parameter too [2l| , one 
finds that the width of the less massive resonance would 
be too small to have escaped experimental detection [26j j . 
Other candidates like the 1650 MeV also have been sug- 
gested [3a |. For the hyperons, we simply chose the split- 
ting between the particles and the respective parity part- 
ners to have the same value as for the nucleons, thus 
keeping the number of parameters to a minimum. This 
assumption agrees quite well with some parity partner 
candidates, like the A(1670) and the £(1750). In the 
case of the 5 the data is unclear. 



III. INCLUSION OF QUARKS 

Following [20| . the effective masses of the quarks are 
generated by the scalar mesons except for a small explicit 
mass term (8m q = 5 MeV and 8m s = 150 MeV for the 
strange quark) and TOq 9 



m* = g q<J a + Sm q + TO 0g , 
m * s = 9s(C + 8ttl s + m 0q , 



(5) 



with g qa = g S Q = 4.0. As was the case for baryons, we 
have also introduced a mass parameter too ? = 200 MeV 
for the quarks. 

We use the Polyakov loop $ as the order parame- 
ter for deconfincment. The field $ is defined via $ = 
■|Tr[exp(i J drA^)], where A4 = iAo is the temporal 
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FIG. 1. (Color online) Mass-radius diagram for star families 
with different temperatures shown for Models A and B (black 
and green lines, respectively). For the cold case we have in- 
cluded a cold crust from [37[ and for the warm case we have 
included a crust with sj ps = 4 in accordance with [38|] . 



component of the SU(3) gauge field [39M4ll | . The cou- 
pling of the quarks to the Polyakov loop is introduced in 
the thermal energy of the quarks (see (2p| for more de- 
tails). All thcrmodynamical quantities are derived from 
the grand canonical potential that includes the effective 
potential [/($,$*, T) , which controls the dynamics of the 
Polyakov-loop. In our approach wc adopt the ansatz pro- 
posed in [42ll43l| 



U 



1 



a(T)$$* 



+ b(T)ln[l - 6$$* + 4($ 3 $* 3 ) - 3($$*) 2 ] , (6) 

with a(T) = a T 4 + a 1 T T 3 + a 2 T$T 2 and b(T) = 6 3 T 3 T. 
The parameters are fixed, as in Ref. j42|, by demanding 
a first order phase transition in the pure gauge sector at 
To = 270 MeV, and that the Stefan-Boltzmann limit of 
a gas of glouns is reached for T — > 00. 

To remove the hadronic contributions from the equa- 
tion of state at high temperatures and densities, we in- 
troduce excluded volume effects. The inclusion of finite- 
volume effects in thermodynamic models for hadronic 
matter w as p roposed in (38l . I44I - I55] . In recent publica- 
tions [20|, [56[ we adopted this ansatz to successfully de- 
scribe a smooth transition from a hadronic to a quark 
dominated system. The modified chemical potential £tj, 
which is connected to the real chemical potential fa of the 
i-th particle species, is obtained by the following relation 



fa = fa - in P , 



(7) 



where = 1 fm 3 is the excluded volume of a particle 
of species i (zero for quarks), and P is the sum over all 
partial pressures. To be thermodynamically consistent, 
all densities (energy density ei, baryon density fa and en- 
tropy density S4) were multiplied by a volume correction 
factor /, which is the ratio of the total volume V and 
the reduced volume V not being occupied. As a con- 
sequence, the chemical potentials of the hadrons are de- 
creased by the quarks, but not vice- versa. In other words, 
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FIG. 2. (Color online) Population (particle density) as a func- 
tion of star radius for a cold star using Model A. Quark den- 
sities are divided by 3. 



as the quarks start appearing, they effectively suppress 
the hadrons by changing their chemical potential, while 
the quarks are only affected through the volume correc- 
tion factor /. 

A shortcoming of the parametrization above is that 
the compressibility of bulk matter at saturation density 
is too high (525 MeV) and the equations of state is a 
bit too stiff a low densities when compared with heavy 
ion collision data from Rcf. 



57] . For this reason, in ad- 



dition to the parametrization discussed above (Model A 
hereafter), we use a second parametrization (Model B) 
that leads to more realistic values for the compressibility 
(k 330 MeV). The values of the modified parameters 
in Model B are v t = 0.5 fm 3 , g Nu = 5.02, gi 1] = -7.79 

and g [ a ] = -0.76 (see Table I of Ref. gfj for the en- 
tire parameter set). The symmetry energy obtained at 
saturation is 32.6 and 32.4 MeV and its slope L is 99 
and 92 MeV, respectively for Model A and B. Both of 
these are in reasonable agreement with symmetry energy 
constraints from Ref. |58[. 

It is important to notice that, although Vi is directed 
related to the position of the phase transitions that oc- 
cur when new particles appear in the system, it cannot 
be varied freely. This parameter is chosen mainly to re- 
produce lattice QCD results at [Ib = for Model A. For 
Model B, one has more freedom for the choice of Vi, but 
the phase transitions reproduced do not change qualita- 
tively with this variation. This happens because in this 
type of unified equation of state, the order in which the 
particles appear depends strongly on charge neutrality. 
For example, negatively charged particles, like the down 
quark, will always be favored with respect to positive 
ones. 



IV. RESULTS 

By assuming chemical equilibrium and charge neutral- 
ity, we obtain an equation of state that describes neu- 
tron star matter. This, in turn, can be used to obtain 



FIG. 3. (Color online) Same as Fig. [2] but for model B. 



a solution of the Tolman-Oppenheimer-Volkoff equation 
[59l l6(i| . which allows us to study how the system be- 
haves under the influence of gravity. We show in Fig. [1] 
sequences of stars with gravitational mass M, and radius 
R for models A and B. For model A, the maximum mass 
obtained for a cold star is 1.96 Mq, which is in agree- 
ment with observations of J1614-2230 (1.97 ± 0.04M Q ) 
[l[ . For a star with temperature of 30 MeV, correspond- 
ing to a proto-neutron star, the maximum mass obtained 
increases to 2.15 Mq. The radius (corresponding to the 
most stable massive star) increases from 10.70 to 15.17 
km, the central baryonic density decreases from 1.16 to 
0.90 fm~ 3 and the central fraction of quarks decreases 
from ~ 60% to ~ 50% when a warm star is considered. 

As for model B, the maximum mass obtained for a cold 
star reduces to 1.64 Mq. For a star with temperature of 
30 MeV the maximum mass obtained is 1.75 Mq. The 
radius of the most stable massive star, in this case, in- 
creases from 9.30 to 14.28 km, the central baryonic den- 
sity decreases from 1.58 to 1.18 fm~ 3 and the central 
fraction of quarks decreases from ~ 50% to ~ 40% when 
a warm star is considered. We clearly see that the price 
paid to obtain a realistic compressibility is the reduction 
of the maximum mass of the neutron star. This effect is 
generally model independent, but is known to be more 
pronounced in the parity model. It is important to keep 
in mind that we have not included in the present calcu- 
lation any rotational or strong magnetic field effects in 
order to keep the effect of the appearance of the parity 
partners more clear. Both rotation and strong magnetic 
fields are known to increase the maximum mass of neu- 
tron stars [6ll, [62j . 

In principle, the correct quantity to be used in order 
to analyze proto-neutron star temporal evolution is the 
entropy per baryon, instead of the temperature. Besides, 
trapped neutrinos should be included in earlier stages 
of star evolution, when the temperature (or entropy per 
baryon) is higher. Although these features are of great 
interest, our goal in this preliminary study of the topic is 
to analyze the application of the model to neutron stars 
in the simplest case of T = and the second simplest case 
of T = 30 MeV. Nevertheless, we chose 30 MeV for the 
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FIG. 4. (Color online) Population (particle density) as a func- 
tion of star radius for a warm star T — 30 MeV using Model 
A. Quark densities are divided by 3. 



finite temperature case because it is approximately the 
maximum temperature reproduced in the center of neu- 
tron stars when the entropy per baryon is fixed to 2. In 
this way, we study the two limiting cases of the star evolu- 
tion (T = and T = 30 MeV). Nevertheless, one can see 
that if the baryon number has to remain constant (like 
the case of isolated stars) not all stars from Fig. [1] exist 
throughout the evolution, and some eventually collapse 
into black wholes. For example, the maximum baryon 
number for parametrization A is 2.73 x 10 57 at T = 
and 2.92 x 10 57 at T = 30, and for parametrization B it 
is 2.24 x 10 57 at T = and 2.32 x 10 57 at T = 30. In a 
detailed study of the matter, all the considerations above 
have to be taken into account. 

The particle population of the cold maximum mass 
neutron stars obtained by both models is shown in Figs. [2] 
and[3l We see that at low densities both stars consist only 
of neutrons, protons, electrons and muons. Towards the 
center of the stars the quarks appear (first the d, followed 
by the u) followed by the chiral partners (of neutrons and 
protons). The threshold for the quarks to appear is 1.56 
po and 2.16 po for parametrizations A and B, respectively. 
It is interesting to notice that in this model the baryons 
appear in small quantities but do not vanish at high den- 
sities. This is a consequence of the suppression of the 
hadrons due to the excluded volume prescription. The 
hadrons are suppressed exponentially by the quark pres- 
sure though not removed entirely. Note, however, that 
their fraction becomes so small that their contribution to 
the thermodynamics of the system is negligible. Further- 
more, at high densities, the amount of nucleons and the 
respective chiral partners is very similar. Particles that 
contain strangeness, despite being included in the model, 
do not appear or appear in very small quantities in the 
cold stars. This is due their large bare masses and the 
slow restoration of the strange scalar (£) field to its zero 
value at high densities. 

The decrease in central density stated above is not 
enough to prevent strange particles to appear in the core 
of warm stars (Figs. 0] and [5]). At low densities we see 
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FIG. 5. (Color online) Same as Fig. [4] but for model B. 



As, £~s, S~s and S°s and at high densities S°s, all with 
their respective chiral partners. The strange quarks only 
appear at densities beyond the ones present in the stud- 
ied stars. Overall, the amount of particles that contain 
strangeness is not high enough to significantly change 
star properties such as mass and radius. These prop- 
erties are influenced mainly by the pressure and energy 
density changing substantially due to the high temper- 
ature, and the Polyakov potential U, which only affects 
the system at finite temperatureQ 

Note that the order of the appearance of the hyperons 
in the star does not happen simply according to their 
vacuum masses. Their appearance is related to the re- 
lation between their effective masses and their effective 
chemical potentials. At zero temperature, for example, a 
particle only appears when its effective chemical poten- 
tial is larger than its effective mass. For finite temper- 
ature, thermal effects also contribute making the above 
requirement more relaxed. But more importantly, the 
order of the particles in a star is strongly constrained by 
charged neutrality. This makes the appearance of nega- 
tive charged hyperons more favorable, followed by charge 
neutral ones. Therefore, the only particles that appear 
in significant quantities for parametrization B at T = 30 
MeV are the E~ followed by the S~ and respective chiral 
partners and they are the ones that still appear at T = 0. 
For parametrization A, because there are no significant 
quantity of hyperons at T = 30 MeV, there are also no 
hyperons in the T = MeV case. 

For the cold case, there are basically two first order 
phase transitions in the star, when the down quarks ap- 
pear and when the up quarks and chiral partners appear. 
For the warm case, both of these transitions are smoother 
and closer to each other. For the symmetric case, rele- 
vant for heavy ion collisions, the phase transitions would 
be stronger for any temperature, since chemical equi- 
librium and charge neutrality push the chiral symmetry 



An example for a Polyakov potential for the deconfinement order 
parameter that also depends on baryon density, i.e. is applicable 
at T = 0, can be found in |(s3l . 
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FIG. 6. (Color online) Thermal evolution of stars of different 
masses for model A. Too denotes the redshifted temperature 
observed at infinity. The observed data consists of circles de- 
noting spin-down ages, squares kinematic ages and diamonds 
the Cas A evolution. 

restoration (related to the appearance of chiral partners) 
to lower densities with respect to the symmetric matter 
case and renders the transition smoother 1261. 



V. THERMAL EVOLUTION 

The thermal evolution of the stars was obtained by 
solving the full set of equations that govern energy bal- 
ance and transport in a relativistic star [64^ . We have 
considered the standard neutrino emission processes for 
both quarks and hadrons, these are the direct Urea (DU), 
modified Urea (MU), and Bremsstrahlung (BM) pro- 
cesses [gH [66]] . The specific heat of the matter is given 
by the traditional specific heat of relativistic fermions, 
as described in Ref. [67j. The thermal conductivity is 
calculated as described in jfi^, |(39[ . Finally, the boundary 
condition that defines the surface temperature of the star 
is discussed in Ref. [7oT - [73j . In our calculations we assume 
a non-magnetized surface, with an accreted envelope of 
AM/A/ = 1.0 X 10~ 9 . 

We have also taken into account neutron superfluid- 
ity when calculating the thermal evolution. The pairing 
patterns considered were the neutron singlet ( 1 <S'o) and 
triplet ( 3 P2) states, as described in [z3. LZ^j . We point 
out that it is p ossible that protons form a singlet super- 
conductor [76lr77j. however, the pairing of protons in the 
core is still not very well understood (see [78| and ref- 
erences therein). For that reason we consider only neu- 
tron superfluidity at this stage. We show in Figs. |5] and 
[7] the cooling of stars with different masses (for Models 
A and B, respectively), being the redshifted surface 
temperature. We also plot a set of prominent observed 
cooling data, where circles denote objects whose age was 
determined through spin-down, squares through the ob- 
ject motion with respect to originating supernova (kine- 
matic a ge) [67], [7!| and diamonds show the evolution of 
Cas A [77(. We note that although our model exhibits 
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FIG. 7. (Color online) Same as Fig. [6] but for model B. 



temperatures and ages comparable with those of Cas A, 
it does not agree with the fast cooling exhibited by this 
object in the last 10 years. We point out, however, that 
differently than in references (76[ we have not removed 
fast neutrino emission processes by assuming that pro- 
tons are in a superconducting state in the entire core. As 
pointed out in [78[ , such stringent condition is more easily 
achieved for low densities stars (less massive), since the- 
oretical models show that proton superconductivity does 
not extend to very high densities. Under these conditions 
(no fast neutrino processes, and for the critical temper- 
ature found in |76|) our model might as well reproduce 
the observed behavior of Cas A. 

The results shown in Fig. [6~] clearly illustrate the effects 
of quark matter and the chiral partners on the cooling of 
the star, which we discuss now. The low density of the 
M/M Q = 1.0 star forbids the presence of chiral partners. 
Consequently the cooling of this object is qualitatively 
the same as that of standard hybrid stars [/J- On the 
other hand, for stars with higher central densities, and 
thus higher masses, we have chiral partner states being 
populated, and the cooling of the star is substantially 
modified. Note that these stars are warmer than their 
low-mass counterparts. Such behavior is the opposite of 
that of ordinary neutron stars, which tend to exhibit a 
faster cooling for heavier objects. In the case of ordinary 
stars the acceleration of the cooling with the increase of 
density is connected to the increasing proton fraction. A 
higher proton fraction throughout the star means a larger 
region where the DU process is present, and thus a faster 
cooling. In our model, however, as the density increases 
(equivalent to smaller radius in Fig. [5]) the proton fraction 
decreases due to the increase of the chiral partners pop- 
ulation. Because of this, in our model, higher mass stars 
present a slower cooling. We also note that the cooling 
of stars in model A and B are qualitatively equal, with 
the difference that Model B cannot reproduce stars with 
masses higher than 1.64 M Q (Fig. [7]). Such an effect of 
slower cooling with increased mass has also been observed 
in hybrid stars with a color-superconducting quark core 
0. 
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FIG. 8. (Color online) Same as in Fig.[6]but for quark matter 
paired with A = 10 MeV. 



We should note here that these results should be con- 
sidered with care. The reason is that so far we have not 
considered neutrino processes involving the chiral part- 
ners, which could certainly affect the cooling. The cross 
section of the processes involving the chiral partners have 
not been determined so far, and are therefore excluded 
from this preliminary study. If one shows that those cross 
sections are much smaller than the ones of the traditional 
particles, i.e. the neutrino luminosity of this process is 
substantially smaller than that of the traditional DU pro- 
cess; then the results shown here would hold. These re- 
sults can as well be easily tested, as soon as there is 
enough statistical data on the relation between pulsar 
age, surface temperature and mass. 

The results of Fig. [6] and [7] show that the cooling of 
hybrid stars in our model is considerably slower than that 
of traditional stars, even the one of the M/M Q = 1.0 star 
that does not contain chiral partners. Additionally, as 
we have already mentioned, we observe that during the 
neutrino dominated era (ages < 10 4 years in ordinary 
objects) more massive stars exhibit a higher temperature. 
This substantially delays the photon cooling era, that 
can be identified by the bending of the cooling tracks at 
later times. We see that our model, with the limitations 
discussed above, agrees fairly well with the cooling data 
of the colder, and older objects, but fails to describe the 
temperature of the very young and hot neutron stars. 

Another possibility to be considered is quark pairing. 
Our study considers a Color-Flavor-Locked (CFL) [81j 
like phase, where all quarks of all colors are paired. Note 
that since we consider small gap values only (A = 10 
MeV ), w hich are not expected to substantially affect the 
EoS [82|, we do not include corrections in the EOS. In 
this way, our analysis is still valid for any quark pair- 
ing scheme (not necessarily color superconductivity), as 
long as it affects all quark flavors in similar way. The 
quarks become a superconductor once the temperature 
falls below the critical temperature T c . For this study we 
assumed, like in reference [7^|, the value of T c = 0.4A. 
Once the quarks become a superconductor the quark DU 
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FIG. 9. (Color online) Same as in Fig.[7]but for quark matter 
paired with A = 10 MeV. 



process is suppressed by a factor e~ A / T , and the modified 
Urea and the Bremsstrahlung process by a factor e~ 2A//T . 
The specific heat of quark matter is also modified by a 
factor 3.2(T C /T) * (2.5 - 1.7(T/T C ) + 3.6(T/T c ) 2 ) e - A / T 

M- 

We stress that the color superconducting phase con- 
sidered here is not, rigorously speaking, the traditional 
CFL phase. In the CFL phase, all quarks of all flavors 
are paired. On another hand, in the 2 flavor supercon- 
ducting state (2SC) one has the BCS pairing between u 
and d quarks of colors r and g, which leaves some cooling 
channels unsuppresscd, thus leading to a faster cooling. 
Differently than the 2SC phase, the model we consider in 
our calculations assumes that u and d quarks of all col- 
ors are paired. As discussed above we consider only small 
gap values (A = 10 MeV), which have little effect to the 
EoS. In that sense our pairing model is similar (but not 
exactly equal) to that discussed in (83[, labeled 2SC+x, 
where we have a 2SC pairing of the u and d of two colors, 
and a residual pairing of the remaining unpaired quark, 
leading to the effective pairing of all u and d quarks. 

We show in Fig. [5] and [§] the cooling of stars whose 
quark cores are in a superconductor state with A = 10 
MeV. One can see that the suppression of the quark emis- 
sion processes allow the stars to be warmer during its 
neutrino cooling era being in better agreement with the 
data of the very young and hot neutron stars. The time 
scale for the onset of the photon cooling era is unchanged. 
As was the case for stars with unpaired quark matter, the 
cooling of stars in models A and B are very similar As 
was the case before, we stress that these results should be 
considered with great care, since the pairing scheme as- 
sumed here is phcnomcnological. A more detailed study, 
taking into account microscopic effects is warranted and 
will be performed in a future study. 
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VI. CONCLUSION 



We applied for the first time the SU(3)f version of the 
parity doublet model to cold, dense, charged neutral and 
chemically equilibrated matter. Note that this approach 
is able to successfully describe the low density as well 
as the high density regime of the QCD phase diagram, 
as it includes a self-consistent deconfinement transition 
to quark matter. With these ingredients, we can study 
for the first time the interplay between the baryon octet, 
their respective chiral partners and quarks in neutron 
stars. The chirally symmetric but still confined matter 
obtained in Ref. |20( disappears when charge neutrality 
and chemical equilibrium are taken into account. This 
is due to the early appearance of the down quark, that 
compensates the positive charge of the proton. Such a 
partial appearance of the quarks was already suggested 
in [84j . The up quark appearance happens, for both cold 
and warm stars, approximately in the same region as the 
chiral partners (signaling the chiral symmetry restora- 
tion) . 

The zero temperature calculation using the same 
paramctrization for the model as the one calibrated for 
low and zero chemical potential (Model A) yields stars 
with masses in agreement with the most massive observed 
pulsar (1.97 ± 0.04M Q This parametrization, how- 
ever, leads to nuclear matter with high compressibility. 
On the other hand, parametrizations leading to more 
realistic compressibility values (like model B) yield less 
massive stars. As already shown in Ref. [26j], corrections 
that account for the baryonic Dirac sea effect such as 
the Relativistic Hartree Approximation (RHA) can im- 
prove this situation. We note that the reason we did not 
consider any extra features that could possible improve 
the situation in this work is because we wanted to study 
in detail the population distribution in the star taking 
into account the relation between baryon chiral partners 
and quarks at high densities, since this had never been 
performed before. Work along this line is in progress. 
The inclusion of finite temperature in the calculation al- 
lows more massive stars in both of the cases studied, but 
still does not qualitatively change the situation presented 
above. 



We have performed cooling simulations for neutron 
stars whose microscopic composition is given by our 
(cold) model. We have found that the presence of the 
chiral partners affects the thermal evolution, effectively 
suppressing the hadronic direct Urea process. This, in 
contrast to other models, yields warmer stars during the 
neutrino cooling era, and delays the onset of the pho- 
ton cooling era. Although we cannot effectively test such 
a prediction at the moment, we will be able to do so, 
hopefully, in the near future. 

We also considered the possibility of pairing in the 
quark phase, where a CFL-like pairing was assumed. The 
suppression of the quark neutrino emission processes, 
brought on by pairing, reduces the total emission of neu- 
trinos further, leading to even warmer stars in the neu- 
trino cooling era. We stress that the cooling results pre- 
sented here should be considered as a first approximation 
only, since there are many factors that still need to be 
considered, like the neutrino emissions from the chiral 
partners, and the effects of stronger quark pairing on the 
EoS, which warrants future work along those lines. We 
believe, however, that the cooling investigation put forth 
in this paper might be a good, qualitative study of the 
cooling of neutron stars within a SU(3)f doublet parity 
model composition, given that the available cooling data 
can be reproduced. Furthermore, a recent work |85( has 
shown that rotation might have an important effect on 
the cooling of compact stars, thus we intend to perform 
2D simulations of the thermal evolution of stars described 
by our model. 
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